[Development of Integrated Database for Experiments and Simulations].

Researchers collect data and use various methods to organize it. Ensuring the reliability and reproducibility of data is crucial, and collaboration across different research fields is on the rise. However, when there is geographical distance, sharing data becomes a challenging task. Therefore, there is a need for the development of a mechanism for sharing data on the web. We have developed an integrated database to facilitate the sharing and management of research data, particularly focusing on small molecules. The integrated database serves as a platform for centralizing data related to small molecules, including their chemical structures, wet lab experimental data, simulation data, and more. It has been constructed as a web application, offering features such as library management for small molecules, registration and viewing of wet lab experiment results, generation of initial conformations for simulations, and data visualization. This enables researchers to efficiently share their research data and collaborate seamlessly, whether within their research group or via cloud-based access that allows project and team members to connect from anywhere. This integrated database plays a critical role in connecting wet lab experiments and simulations, enabling researchers to cross-reference and analyze experimental data comprehensively. It serves as an essential tool to advance research and foster idea generation.

[Development of Membrane Permeability Coefficient by Means of Novel Molecular Dynamics Methods].

The membrane permeability, and its evaluation, is crucial factor in the process of uptake of compounds from outside to inside the cell and in the inhibition of the activity of disease-causing target proteins. Although molecular dynamics (MD) simulations have been shown to be able to reproduce the conformational changes of compounds occurring during membrane permeation, it is still challenging to extract the membrane permeability at an affordable computational workload solely by conventional MD. Indeed, the time scale accessible by MD is far below the one characterizing the actual permeation process. Phenomena occurring in living organisms escaping the reach of standard MD are generally referred to as biological rare events, and the membrane permeation process is one of them. To overcome this time-scale problem, several enhanced sampling methods have been proposed over the years to improve conformational sampling. In this review, a hybrid sampling method that combines the parallel cascade selection MD (PaCS-MD) and the outlier flooding method (OFLOOD), introduced and developed by our group, is proposed as a tool to study the membrane permeation from structural sampling (rare-event sampling). The obtained trajectories are used to estimate the free energy profiles for the membrane permeation and to compute the membrane permeation coefficients. Moreover, we present an example of application of the free energy reaction network method as a versatile way for incorporating explicitly into reaction coordinates the degrees of freedom related to internal motion.

Development of an FeII Complex Exhibiting Intermolecular Proton Shifting Coupled Spin Transition.

In this study, we investigated reversible intermolecular proton shifting (IPS) coupled with spin transition (ST) in a novel FeII complex. The host FeII complex and the guest carboxylic acid anion were connected by intermolecular hydrogen bonds (IHBs). We extended the intramolecular proton transfer coupled ST phenomenon to the intermolecular system. The dynamic phenomenon was confirmed by variable-temperature single-crystal X-ray diffraction, neutron crystallography, and infrared spectroscopy. The mechanism of IPS was further validated using density functional theory calculations. The discovery of IPS-coupled ST in crystalline molecular materials provides good insights into fundamental processes and promotes the design of novel multifunctional materials with tunable properties for various applications, such as optoelectronics, information storage, and molecular devices.

Grinding-Induced Water Solubility Exhibited by Mechanochromic Luminescent Supramolecular Fibers.

Most mechanochromic luminescent compounds are crystalline and highly hydrophobic; however, mechanochromic luminescent molecular assemblies comprising amphiphilic molecules have rarely been explored. This study investigated mechanochromic luminescent supramolecular fibers composed of dumbbell-shaped 9,10-bis(phenylethynyl)anthracene-based amphiphiles without any tetraethylene glycol (TEG) substituents or with two TEG substituents. Both amphiphiles formed water-insoluble supramolecular fibers via linear hydrogen bond formation. Both compounds acquired water solubility when solid samples composed of supramolecular fibers are ground. Grinding induces the conversion of 1D supramolecular fibers into micellar assemblies where fluorophores can form excimers, thereby resulting in a large redshift in the fluorescence spectra. Excimer emission from the ground amphiphile without TEG chains is retained after dissolution in water. The micelles are stable in water because hydrophilic dendrons surround the hydrophobic luminophores. By contrast, when water is added to a ground amphiphile having TEG substituents, fragmented supramolecular fibers with the same molecular arrangement as the initial supramolecular fibers are observed, because fragmented fibers are thermodynamically preferable to micelles as the hydrophobic arrays of fluorophores are covered with hydrophilic TEG chains. This leads to the recovery of the initial fluorescent properties for the latter amphiphile. These supramolecular fibers can be used as practical mechanosensors to detect forces at the mesoscale.

FMO-guided design of darunavir analogs as HIV-1 protease inhibitors.

The prevalence of HIV-1 infection continues to pose a significant global public health issue, highlighting the need for antiretroviral drugs that target viral proteins to reduce viral replication. One such target is HIV-1 protease (PR), responsible for cleaving viral polyproteins, leading to the maturation of viral proteins. While darunavir (DRV) is a potent HIV-1 PR inhibitor, drug resistance can arise due to mutations in HIV-1 PR. To address this issue, we developed a novel approach using the fragment molecular orbital (FMO) method and structure-based drug design to create DRV analogs. Using combinatorial programming, we generated novel analogs freely accessible via an on-the-cloud mode implemented in Google Colab, Combined Analog generator Tool (CAT). The designed analogs underwent cascade screening through molecular docking with HIV-1 PR wild-type and major mutations at the active site. Molecular dynamics (MD) simulations confirmed the assess ligand binding and susceptibility of screened designed analogs. Our findings indicate that the three designed analogs guided by FMO, 19-0-14-3, 19-8-10-0, and 19-8-14-3, are superior to DRV and have the potential to serve as efficient PR inhibitors. These findings demonstrate the effectiveness of our approach and its potential to be used in further studies for developing new antiretroviral drugs.

Manipulation of photosynthetic energy transfer by vibrational strong coupling.

Uncovering the mystery of efficient and directional energy transfer in photosynthetic organisms remains a critical challenge in quantum biology. Recent experimental evidence and quantum theory developments indicate the significance of quantum features of molecular vibrations in assisting photosynthetic energy transfer, which provides the possibility of manipulating the process by controlling molecular vibrations. Here, we propose and theoretically demonstrate efficient manipulation of photosynthetic energy transfer by using vibrational strong coupling between the vibrational state of a Fenna-Matthews-Olson (FMO) complex and the vacuum state of an optical cavity. Specifically, based on a full-quantum analytical model to describe the strong coupling effect between the optical cavity and molecular vibration, we realize efficient manipulation of energy transfer efficiency (from 58% to 92%) and energy transfer time (from 20 to 500 ps) in one branch of FMO complex by actively controlling the coupling strength and the quality factor of the optical cavity under both near-resonant and off-resonant conditions, respectively. Our work provides a practical scenario to manipulate photosynthetic energy transfer by externally interfering molecular vibrations via an optical cavity and a comprehensible conceptual framework for researching other similar systems.

Machine learning-guided design of potent darunavir analogs targeting HIV-1 proteases: A computational approach for antiretroviral drug discovery.

In the pursuit of novel antiretroviral therapies for human immunodeficiency virus type-1 (HIV-1) proteases (PRs), recent improvements in drug discovery have embraced machine learning (ML) techniques to guide the design process. This study employs ensemble learning models to identify crucial substructures as significant features for drug development. Using molecular docking techniques, a collection of 160 darunavir (DRV) analogs was designed based on these key substructures and subsequently screened using molecular docking techniques. Chemical structures with high fitness scores were selected, combined, and one-dimensional (1D) screening based on beyond Lipinski's rule of five (bRo5) and ADME (absorption, distribution, metabolism, and excretion) prediction implemented in the Combined Analog generator Tool (CAT) program. A total of 473 screened analogs were subjected to docking analysis through convolutional neural networks scoring function against both the wild-type (WT) and 12 major mutated PRs. DRV analogs with negative changes in binding free energy ( ΔΔ G bind ) compared to DRV could be categorized into four attractive groups based on their interactions with the majority of vital PRs. The analysis of interaction profiles revealed that potent designed analogs, targeting both WT and mutant PRs, exhibited interactions with common key amino acid residues. This observation further confirms that the ML model-guided approach effectively identified the substructures that play a crucial role in potent analogs. It is expected to function as a powerful computational tool, offering valuable guidance in the identification of chemical substructures for synthesis and subsequent experimental testing.

Latrunculin resistance mechanism of non-conventional actin NAP1 uncovered by molecular dynamics simulations.

Monomeric G-actin polymerizes into F-actin to perform various cellular functions. Actin depolymerization drugs, such as latrunculin-A (Lat-A), inhibit filament formation and disrupt the cytoskeleton. Interestingly, the green algae Chlamydomonas alternatively produces a non-conventional actin, NAP1, that responds to inhibition by latrunculin. However, the molecular mechanism underlying latrunculin resistance of NAP1 remains unclear because of the difficulty due to its low in vitro polymerizability. Instead of biochemical experiments, we performed molecular dynamics (MD) simulations to investigate whether NAP1 has a lower affinity for Lat-A than the conventional actins. Our phylogenetic comparison of the binding free energies shows that Lat-A is evolutionarily optimized for skeletal muscles. By decomposing the binding free energy into each amino acid residue, we found that some residues in NAP1 play an important role in latrunculin resistance, suggesting that the primary mechanism of latrunculin resistance is the loss of affinity for Lat-A due to substitutions. In conclusion, our binding-free-energy calculations using MD simulations provide the critical insight that loss of affinity is the direct mechanism of latrunculin resistance.

Synthesis and biological evaluation of 2'-hydroxychalcone derivatives as AMPK activators.

A series of 2'-hydroxychalcone derivatives with various substituents on B-ring were synthesized and evaluated for AMP-activated protein kinase (AMPK) activation activity in podocyte cells. The results displayed that hydroxy, methoxy and methylenedioxy groups on B-ring could enhance the activitiy better than O-saturated alkyl, O-unsaturated alkyl or other alkoxy groups. Compounds 27 and 29 possess the highest fold change of 2.48 and 2.73, respectively, which were higher than those of reference compound (8) (1.28) and metformin (1.88). Compounds 27 and 29 were then subjected to a concentration-response study to obtain the EC50 values of 2.0 and 4.8 µM, respectively and MTT assays also showed that cell viability was not influenced by the exposure of podocytes to compounds 27 and 29 at concentrations up to 50 µM. In addition, compound 27 was proved to activate AMPK via calcium/calmodulin-dependent protein kinase kinase ß (CaMKKß)-dependent pathway without affecting intracellular calcium levels. The computational study showed that the potent compounds exhibited stronger ligand-binding strength to CaMKKß, particularly compounds 27 (-8.4 kcal/mol) and 29 (-8.0 kcal/mol), compared to compound 8 (-7.5 kcal/mol). Fragment molecular orbital (FMO) calculation demonstrated that compound 27 was superior to compound 29 due to the presence of methyl group, which amplified the binding by hydrophobic interactions. Therefore, compound 27 would represent a promising AMPK activator for further investigation of the treatment of diabetes and diabetic nephropathy.

Pathological mutations promote proteolysis of mitochondrial tRNA-specific 2-thiouridylase 1 (MTU1) via mitochondrial caseinolytic peptidase (CLPP).

MTU1 controls intramitochondrial protein synthesis by catalyzing the 2-thiouridine modification of mitochondrial transfer RNAs (mt-tRNAs). Missense mutations in the MTU1 gene are associated with life-threatening reversible infantile hepatic failure. However, the molecular pathogenesis is not well understood. Here, we investigated 17 mutations associated with this disease, and our results showed that most disease-related mutations are partial loss-of-function mutations, with three mutations being particularly severe. Mutant MTU1 is rapidly degraded by mitochondrial caseinolytic peptidase (CLPP) through a direct interaction with its chaperone protein CLPX. Notably, knockdown of CLPP significantly increased mutant MTU1 protein expression and mt-tRNA 2-thiolation, suggesting that accelerated proteolysis of mutant MTU1 plays a role in disease pathogenesis. In addition, molecular dynamics simulations demonstrated that disease-associated mutations may lead to abnormal intermolecular interactions, thereby impairing MTU1 enzyme activity. Finally, clinical data analysis underscores a significant correlation between patient prognosis and residual 2-thiolation levels, which is partially consistent with the AlphaMissense predictions. These findings provide a comprehensive understanding of MTU1-related diseases, offering prospects for modification-based diagnostics and novel therapeutic strategies centered on targeting CLPP.

Structural and thermodynamic insights into antibody light chain tetramer formation through 3D domain swapping.

Overexpression of antibody light chains in small plasma cell clones can lead to misfolding and aggregation. On the other hand, the formation of amyloid fibrils from antibody light chains is related to amyloidosis. Although aggregation of antibody light chain is an important issue, atomic-level structural examinations of antibody light chain aggregates are sparse. In this study, we present an antibody light chain that maintains an equilibrium between its monomeric and tetrameric states. According to data from X-ray crystallography, thermodynamic and kinetic measurements, as well as theoretical studies, this antibody light chain engages in 3D domain swapping within its variable region. Here, a pair of domain-swapped dimers creates a tetramer through hydrophobic interactions, facilitating the revelation of the domain-swapped structure. The negative cotton effect linked to the ß-sheet structure, observed around 215 nm in the circular dichroism (CD) spectrum of the tetrameric variable region, is more pronounced than that of the monomer. This suggests that the monomer contains less ß-sheet structures and exhibits greater flexibility than the tetramer in solution. These findings not only clarify the domain-swapped structure of the antibody light chain but also contribute to controlling antibody quality and advancing the development of future molecular recognition agents and drugs.

A sulfonamide chalcone inhibited dengue virus with a potential target at the SAM-binding site of viral methyltransferase.

Dengue infection is a global health problem as climate change facilitates the spread of mosquito vectors. Infected patients could progress to severe plasma leakage and hemorrhagic shock, where current standard treatment remains supportive. Previous reports suggested that several flavonoid derivatives inhibited mosquito-borne flaviviruses. This work aimed to explore sulfonamide chalcone derivatives as dengue inhibitors and to identify molecular targets. We initially screened 27 sulfonamide chalcones using cell-based antiviral and cytotoxic screenings. Two potential compounds, SC22 and SC27, were identified with DENV1-4 EC50s in the range of 0.71-0.94 and 3.15-4.46 µM, and CC50s at 14.63 and 31.02 µM, respectively. The compounds did not show any elevation in ALT or Cr in C57BL/6 mice on the 1st, 3rd, and 7th days after being administered intraperitoneally with 50 mg/kg SC22 or SC27 in a single dose. Moreover, the SAM-binding site of NS5 methyltransferase was a potential target of SC27 identified by computational and enzyme-based assays. The main target of SC22 was in a late stage of viral replication, but the exact target molecule had yet to be identified. In summary, a sulfonamide chalcone, SC27, was a potential DENV inhibitor that targeted viral methyltransferase. Further investigation should be the study of the structure-activity relationship of SC27 derivatives for higher potency and lower toxicity.

Design of electron-donating group substituted 2-PAM analogs as antidotes for organophosphate insecticide poisoning.

The use of organophosphate (OPs) pesticides is widespread in agriculture and horticulture, but these chemicals can be lethal to humans, causing fatalities and deaths each year. The inhibition of acetylcholinesterase (AChE) by OPs leads to the overstimulation of cholinergic receptors, ultimately resulting in respiratory arrest, seizures, and death. Although 2-pralidoxime (2-PAM) is the FDA-approved drug for treating OP poisoning, there is difficulty in blood-brain barrier permeation. To address this issue, we designed and evaluated a series of 2-PAM analogs by substituting electron-donating groups on the para and/or ortho positions of the pyridinium core using in silico techniques. Our PCM-ONIOM2 (MP2/6-31G*:PM7//B3LYP/6-31G*:UFF) binding energy results demonstrated that 13 compounds exhibited higher binding energy than 2-PAM. The analog with phenyl and methyl groups substituted on the para and ortho positions, respectively, showed the most favorable binding characteristics, with aromatic residues in the active site (Y124, W286, F297, W338, and Y341) and the catalytic residue S203 covalently bonding with paraoxon. The results of DS-MD simulation revealed a highly favorable apical conformation of the potent analog, which has the potential to enhance reactivation of AChE. Importantly, newly designed compound demonstrated appropriate drug-likeness properties and blood-brain barrier penetration. These results provide a rational guide for developing new antidotes to treat organophosphate insecticide toxicity.

Structure-yeast α-glucosidase inhibitory activity relationship of 9-O-berberrubine carboxylates.

Thirty-five 9-O-berberrubine carboxylate derivatives were synthesized and evaluated for yeast α-glucosidase inhibitory activity. All compounds demonstrated better inhibitory activities than the parent compounds berberine (BBR) and berberrubine (BBRB), and a positive control, acarbose. The structure-activity correlation study indicated that most of the substituents on the benzoate moiety such as methoxy, hydroxy, methylenedioxy, benzyloxy, halogen, trifluoromethyl, nitro and alkyl can contribute to the activities except multi-methoxy, fluoro and cyano. In addition, replacing benzoate with naphthoate, cinnamate, piperate or diphenylacetate also led to an increase in inhibitory activities except with phenyl acetate. 9, 26, 27, 28 and 33 exhibited the most potent α-glucosidase inhibitory activities with the IC50 values in the range of 1.61-2.67 µM. Kinetic study revealed that 9, 26, 28 and 33 interacted with the enzyme via competitive mode. These four compounds were also proved to be not cytotoxic at their IC50 values. The competitive inhibition mechanism of these four compounds against yeast α-glucosidase was investigated using molecular docking and molecular dynamics simulations. The binding free energy calculations suggest that 26 exhibited the strongest binding affinity, and its binding stability is supported by hydrophobic interactions with D68, F157, F158 and F177. Therefore, 9, 26, 28 and 33 would be promising candidates for further studies of antidiabetic activity.

Deciphering the Potential of Multidimensional Carbon Materials for Surface-Enhanced Raman Spectroscopy through Density Functional Theory.

Surface-enhanced Raman spectroscopy (SERS) is a potent analytical tool, particularly for molecular identification and structural analysis. Conventional metallic SERS substrates, however, suffer from low reproducibility and compatibility with biological molecules. Recently, metal-free SERS substrates based on chemical enhancement have emerged as a promising alternative with carbon-based materials offering excellent reproducibility and compatibility. Nevertheless, our understanding of carbon materials in SERS remains limited, which hinders their rational design. Here we systematically explore multidimensional carbon materials, including zero-dimensional fullerenes (C60), one-dimensional carbon nanotubes, two-dimensional graphene, and their B-, N-, and O-doped derivatives, for SERS applications. Using density functional theory, we elucidate the nonresonant polarizability-enhanced and resonant charge-transfer-based chemical enhancement mechanisms of these materials by evaluating their static/dynamic polarizability and electron excitation properties. This work provides a critical reference for the future design of carbon-based SERS substrates, opening a new avenue in this field.

The 8-bromobaicalein alleviated chikungunya-induced musculoskeletal inflammation and reduced the viral load in healthy adult mice.

Chikungunya virus is a re-emerging arbovirus that has caused epidemic outbreaks in recent decades. Patients in older age groups with high viral load and severe immunologic response during acute infection are likely to develop chronic arthritis and severe joint pain. Currently, no antiviral drug is available. Previous studies suggested that a flavone derivative, 8-bromobaicalein, was a potential dengue and Zika replication inhibitor in a cell-based system targeting flaviviral polymerase. Here we characterized that 8-bromobaicalein inhibited chikungunya virus replication with EC50 of 0.49 ± 0.11 µM in Vero cells. The molecular target predicted at viral nsP1 methyltransferase using molecular binding and fragment molecular orbital calculation. Additionally, oral administration of 250 mg/kg twice daily treatment alleviated chikungunya-induced musculoskeletal inflammation and reduced viral load in healthy adult mice. Pharmacokinetic analysis indicated that the 250 mg/kg administration maintained the compound level above EC99.9 for 12 h. Therefore, 8-bromobaicalein should be a potential candidate for further development as a pan-arboviral drug.

Origin of Homochirality in Amino Acids Induced by Lyman-α Irradiation in the Early Stage of the Milky Way.

The enantiomeric excess (ee) of l-form amino acids found in the Murchison meteorite poses some issues about the cosmic origin of their chirality. Circular dichroism (CD) spectra of amino acids in the far-ultraviolet (FUV) at around 6.8 eV (182 nm) indicate that the circularly polarized light can induce ee through photochemical reactions. Here, we resort to ab initio calculations to extract the CD spectra up to the vacuum-ultraviolet (VUV) region (∼11 eV), and we propose a novel equation to compute the ee applicable to a wider range of light frequency than what is available to date. This allows us to show that the strength of the induced ee (|ee|) in the 10 eV VUV region is comparable to the one in the 6.8 eV FUV region. This feature is common for some key amino acids (alanine, 2-aminobutyric acid, and valine). In space, intense Lyman-α (Lyα) light of 10.2 eV is emitted from star forming regions. This study provides a theoretical basis that Lyα emitter from an early starburst in the Milky Way plays a crucial role in initiating the ee of amino acids.

New insights into the oxidation process from neutron and X-ray crystal structures of an O2-sensitive [NiFe]-hydrogenase.

[NiFe]-hydrogenase from Desulfovibrio vulgaris Miyazaki F is an O2-sensitive enzyme that is inactivated in the presence of O2 but the oxidized enzyme can recover its catalytic activity by reacting with H2 under anaerobic conditions. Here, we report the first neutron structure of [NiFe]-hydrogenase in its oxidized state, determined at a resolution of 2.20 Å. This resolution allowed us to reinvestigate the structure of the oxidized active site and to observe the positions of protons in several short hydrogen bonds. X-ray anomalous scattering data revealed that a part of the Ni ion is dissociated from the active site Ni-Fe complex and forms a new square-planar Ni complex, accompanied by rearrangement of the coordinated thiolate ligands. One of the thiolate Sγ atoms is oxidized to a sulfenate anion but remains attached to the Ni ion, which was evaluated by quantum chemical calculations. These results suggest that the square-planar complex can be generated by the attack of reactive oxygen species derived from O2, as distinct from one-electron oxidation leading to a conventional oxidized form of the Ni-Fe complex. Another major finding of this neutron structure analysis is that the Cys17S thiolate Sγ atom coordinating to the proximal Fe-S cluster forms an unusual hydrogen bond with the main-chain amide N atom of Gly19S with a distance of 3.25 Å, where the amide proton appears to be delocalized between the donor and acceptor atoms. This observation provides insight into the contribution of the coordinated thiolate ligands to the redox reaction of the Fe-S cluster.

Development of an Iron(II) Complex Exhibiting Thermal- and Photoinduced Double Proton-Transfer-Coupled Spin Transition in a Short Hydrogen Bond.

Multiple proton transfer (PT) controllable by external stimuli plays a crucial role in fundamental chemistry, biological activity, and material science. However, in crystalline systems, controlling multiple PT, which results in a distinct protonation state, remains challenging. In this study, we developed a novel tridentate ligand and iron(II) complex with a short hydrogen bond (HB) that exhibits a PT-coupled spin transition (PCST). Single-crystal X-ray and neutron diffraction measurements revealed that the positions of the two protons in the complex can be controlled by temperature and photoirradiation based on the thermal- and photoinduced PCST. The obtained results suggest that designing molecules that form short HBs is a promising approach for developing multiple PT systems in crystals.

Identification of Promising Sulfonamide Chalcones as Inhibitors of SARS-CoV-2 3CLpro through Structure-Based Virtual Screening and Experimental Approaches.

3CLpro is a viable target for developing antiviral therapies against the coronavirus. With the urgent need to find new possible inhibitors, a structure-based virtual screening approach was developed. This study recognized 75 pharmacologically bioactive compounds from our in-house library of 1052 natural product-based compounds that satisfied drug-likeness criteria and exhibited good bioavailability and membrane permeability. Among these compounds, three promising sulfonamide chalcones were identified by combined theoretical and experimental approaches, with SWC423 being the most suitable representative compound due to its competitive inhibition and low cytotoxicity in Vero E6 cells (EC50 = 0.89 ± 0.32 µM; CC50 = 25.54 ± 1.38 µM; SI = 28.70). The binding and stability of SWC423 in the 3CLpro active site were investigated through all-atom molecular dynamics simulation and fragment molecular orbital calculation, indicating its potential as a 3CLpro inhibitor for further SARS-CoV-2 therapeutic research. These findings suggested that inhibiting 3CLpro with a sulfonamide chalcone such as SWC423 may pave the effective way for developing COVID-19 treatments.